Localization transition on complex networks via spectral statistics 
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The spectral statistics of complex networks are numerically studied. The features of the Anderson 
metal-insulator transition are found to be similar for a wide range of different networks. A metal- 
insulator transition as a function of the disorder can be observed for different classes of complex 
networks for which the average connectivity is small. The critical index of the transition corresponds 
to the mean field expectation. When the connectivity is higher, the amount of disorder needed to 
reach a certain degree of localization is proportional to the average connectivity, though a precise 
transition cannot be identified. The absence of a clear transition at high connectivity is probably 
due to the very compact structure of the highly connected networks, resulting in a small diameter 
even for a large number of sites. 

PACS numbers: 89.75.Hc,72.15.Rn 



The Anderson transition predicts a transition from ex- 
tended (metallic) to localized (insulating) eigenstates as 
a function of the disorder or energy of a quantum sys- 
tem Ij. This second-order phase transition turns out to 
be a very general property related to the transport of 
quantum and classical waves in disordered systems and 
signatures of it have been observed for electrons in met- 
als, microwaves in waveguides, light in liquids and gels 
and acoustical waves in the earth crust 2]. For all these 
systems the clearest signature of the transition is the dif- 
ferent spread of a signal injected into the system. While 
in a metallic phase the injected wave will spread all over 
the system, in the insulating phase it will be localized in 
the vicinity of the injection point. This effect is a result of 
the constructive interference between time-reversed path 
throughout the system. Thus, the details of the tran- 
sition are strongly influenced by the dimensionality and 
topology of the system. 

The lower critical dimension, below which the system 
is localized for all values of disorder, is believed to be two 
[sl, since the probability of returning to the origin (i.e., 
constructive interference due to time reversal symmetr 
is finite below d = 2. The upper critical dimension 
(beyond which the critical exponents reach their mean- 
field values) remains uncertain although it is argued to 
be infinity 5] . The parameters defining the transition are 
traditionally given as the critical disorder Wc expressed 
in terms of the width of the distribution from which the 
on-sites energies in the Anderson model are drawn, and 
the critical exponent ly (for definitions see Sec. EJ. For 
square lattices with dimensionality d = 3, 4 the values of 
Wc and v are well established : for d = 3, Wc ~ 16.5 and 

~ 1-5 00, while for d = 4, ~ 35 and - 1 
For higher dimensions the following extrapolation was of- 
fered 9]: ly - 0.8/(d -2) + 0.5 and Wc - 16.5(fi - 2), 
which was obtained by studying the transition on bifrac- 
tal topologies. Thus, the mean field critical index value 
of 1/2 is obtained in the the upper critical dimension 
d = oo. 

There has been recently much interest in the proper- 
ties of random scale-free networks 0, 0, [T3. fl^ 



ITsl IT^ IT^ ITM IT^ . These networks are characterized 
by the fact that each node is connected with some finite 
probability to any other node in the graph, which is very 
different from the usual topology of a real space lattice in 
which nodes are connected only to their neighbors. This 
leads to a very interesting behavior of the graphs when 
properties such as percolation, cluster structures, paths 
length etc. are considered ^Jl- Interest in the influence 
this unusual topology has on the properties of wave inter- 
ference (i.e., the Anderson transition of these graphs) is 
rising. Indeed, recently the Anderson transition in partic- 
ular networks, namely the small- world networks pol l2ll | 
and the Caylcy tree ^21 were studied. Older work on the 
localization properties of sparse random matrices is 
directly relevant to Erdos-Renyi graphs. Since scale free 
networks have an unusual topology for which anomalous 
classical properties have been found 0,|2^ it is of partic- 
ular interest to study the Anderson localization in these 
networks. 

Essentially, the probability to return to the origin de- 
fines the dimensionality of a system for the Anderson 
transition. Therefore, one may speculate that random 
graphs, which have only very long closed trajectories, 
correspond to systems with an infinite dimension |2(i| |. 
On the other hand, the critical disorder which depends 
roughly on the number of nearest neighbors Z is expected 
to follow Q Wc Z , which for a random graph with an 
average degree (k) (i.e., the average number of connec- 
tions per node) corresponds to Wc (k). Thus, one may 
expect here an interesting situation in which the criti- 
cal index ly, which is determined by the dimensionality is 
close to a half, while the critical disorder is determined by 
(k) . This is very different than the situation described by 
the extrapolation given above, where the critical disorder 
for an infinite dimension should also be infinite. 

Beyond the general interest in investigating the Ander- 
son transition on scale free networks, and the fresh out- 
look it might provide on the localization phenomenon, 
the metal-insulator transition can provide insights into 
the functionality of complex networks. Consider for ex- 
ample an optical communication network. In such a net- 
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work the edges of the graph represent optical fibers in 
which light propagates and the nodes represents a beam 
splitter which redistributes the incoming wave into the 
outgoing bonds connected to the node. Since for high 
quality optical fiber there are essentially no losses or de- 
coherence on the bonds, the amplitude of the transmitted 
wave will not depend on the bond length (on the other 
hand phase will depend on the length). This network 
may be mapped on a tight binding Hamiltonian of the 
type described in Sec. El ^3- interesting question for 
such a scale free network is whether a wave injected into 
one of the nodes will produce a signal at all other node 
(which in the language of the Anderson localization is 
equivalent to the question is the system metallic) or only 
at a finite set of other nodes (i.e., the system is insulat- 
ing). A metal-insulator transition in this network corre- 
sponds to a phase transition between those two phases 
as function of the properties of the nodes. Generally, for 
any complex network in which information propagates in 
a wave-like fashion and interference is possible, the An- 
derson transition will limit the spread of the information 
throughout the network. 

The paper is organized in the following way: In the 
next section (Sec. we describe the different networks 
which were considered, while in Sec|n]the spectral statis- 
tics method which was applied in order to identify the 
metal-insulator transition is outlined. In Sec. IIIII the re- 
sults are depicted and some general characteristics of the 
localization on complex networks are discussed. 



I. CHARACTERISTICS OF THE DIFFERENT 
NETWORKS 

Our main goal is to study the Anderson transition for 
different complex networks. In this section we shall define 
the characteristics of the networks which will be consid- 
ered. 



follows the Poisson distribution: 

fc! ■ 



P(fc) 



(1) 



The diameter of such a graph follows: d ^ \nN, 
similar to a random graph. In Sec. IIIII we have 
specifically calculated the level distribution for (fc) ~ 
3,3.1,3.2,3.5,4,5,7.5 and 10. 



C. Scale-Free Networks 

Scale free (SF) networks T^l are networks where the 
degree distribution (i.e., fraction of sites with k connec- 
tions) decays as a power-law. The degree distribution is 



given by 



cays i 

m ■■ 



p{k) 



ck < k < K, 



where c = (A - l)m^^^ and K = miV^/(^~^' 14], A is 
the power-law exponent, to is a lower cutoff, and K is 
an upper cutoff . Thus, there are no sites with degree 
below TO or above K. The diameter of the SF network 
can be regarded as the mean distance of the sites from 
the site with the highest degree. For graphs with 2 < 
A < 3 the distance behaves as d ~ ln(ln(A^)) J^], and 
for A = 3 as d - In(iV)/ ln(ln(A^)) This anomalous 
behavior stems from the structure of the network where a 
small core containing most of the high degree sites has a 
very small diameter. For higher values of A the distance 
behaves as in ER, i.e., d ~ In(iV). The (fc) of a SF graph 
is obtained by the following expression: 



(fc) 



1-A K 



2-A 



,2-A 



2-A K^-^-m 



,1-A 



(2) 



A. Random Graph 



For A > 2 and large enough N the average degree ,(fc), 
is a constant. The SF networks analyzed in this paper 
correspond to A = 3.5,4,5 with to = 2 (lower cutoff), 
and A = 4, 5 with to = 3. Due to their small diameter, 
SF networks with A < 3 were not considered. 



A random graph (or- random regular graph) is a graph 
with N nodes, each is connected to exactly fc random 
neighbors ^Jl ■ The diameter of a graph is the maximal 
distance between any pair of its nodes. In a random 
graph the diameter d is proportional to IniV. In Sec. IIIII 
we shall present results of the level spacing distribution 
for random-regular graphs with fc = 3 . 



B. Erdos-Renyi Graphs 

In their classical model from 1959 Erdos and Renyi 
(ER) 28] describe a graph with N nodes where every pair 
of nodes is connected with probability p resulting in (fc) ~ 
Np. For a large random graph the degree distribution 



D. Double peaked distributions 

In order to find hierarchical relation between the dif- 
ferent graphs, we studied also some variations on these 
graphs. For a random graph we changed the degree of 
a small percent of the nodes, so we have a graph with 
double peaked distribution. Thus, the average connec- 
tivity, (fc), is the average degree of the nodes. Several 
examples were taken: changing 5% of the nodes to fc = 5 
(instead of fc = 3) resulting in (fc) = 3.1, or changing 5% 
of the nodes to fc = 10 ((fc) = 3.35). Replacing 20% of 
the nodes connectivity for the previous cases will result 
in (fc) = 3.4 (for fc = 5 nodes) and (fc) = 4.4 (for fc = 10 
nodes). Additionally, in order to relate with previous re- 
sults of the metal insulator transition on a Cayley-tree 
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[23 | we checked a tree in which 5% of its nodes have 
higher degree (fc = 4) resulting in an average connectiv- 
ity 3.05 and creating few closed trajectories - loops. 



II. METHOD 

Now we turn to the calculation of the spectral statis- 
tics of these networks. First, one must construct the 
appropriate network structure, i.e., to determine which 
node is connected to which. This is achieved using the 
following algorithm 0, : 

1. For each site choose a degree from the required 
distribution. 

2. Create a list in which each site is repeated as many 
times as its degree. 

3. Choose randomly two sites from the list and connect 
this pair of site as long as they are different sites. 

4. Remove the pair from the list. Return to 3. 

The diameter of a graph is calculated by building shells 
of sites [23. The inner shell contains the node with the 
highest degree, the next contains all of its neighbors, and 
so on. Of course, each node is counted only once. The 
diameter of the system is then determined by the number 
of shells. Two more options which were considered are 
defining the diameter by the most highly populated shell, 
or by averaging over the shells. The diameter obtained 
by the various methods are quite similar. 

The energy spectrum is calculated using the usual 
tight-binding Hamiltonian, 



H 



(3) 



where first term of the Hamiltonian stands for the dis- 
ordered on-site potential on each node i of the network. 
The on-site energies, are uniformly distributed over the 
range —W/2 < < 14^/2. The second term corresponds 
to the hopping matrix element which is set to 1, and 
denotes nearest neighbor nodes which are determined ac- 
cording to the network structure. We diagonalize the 
Hamiltonian exactly, and obtain N eigenvalues Ei (where 
iV is the number of nodes in the graph) and eigenvectors 
ipi. Then we calculate the distribution P{s) of adjacent 
level spacings s, where s = (-Ei+i — Ei)/{Ei+i — Ei), and 
(. . .) denotes averaging over different realizations of dis- 
order and when relevant also over different realizations 
of node connectivities. 

One expects the distribution to shift as function of 
the on-site disorder from the Wigner surmise distribution 
(characteristic of extended states) : 



Pwis) 



ITS 



■ exp 



TTS 



(4) 



at weak disorder to a Poisson distribution (characteristic 
of localized states) at strong disorder: 




FIG. 1: The distribution P(s) for a 500 sites scale-free graph 
with A = 4 (m = 2). A clear transition from Wigner to 
Poisson is observed as a function of disorder. 



An example for such a transition is presented in Fig. ^ 
where a scale-free graph with A = 4 and m — 2 was con- 
sidered. As W increases P{s) shifts toward the Poisson 
distribution. Additional hallmark features of the Ander- 
son transition, such as the fact that all curves intersects 
at s = 2 and the peak of the distribution " climbs" along 
the Poisson curve for larger values of W are also appar- 
ent. Similar transition from Wigner to Poisson statis- 
tics is seen also for the other networks considered in this 
study. 

The transition point can be determined more accu- 
rately from calculating _30j 
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/2°° P{s)ds - /2°° Pw{s)ds 
Pp{s)ds - Pio{s)ds ' 



(6) 



where 7 as the distribution tends toward the Wigner 
distribution, and 7 ^ 1 if the distribution approaches the 
Poisson distribution. One expects that as the system size 
increases, the finite size corrections will become smaller 
resulting in a distribution closer to a Wigner distribution 
in the metallic regime and to Poisson in the localized 
one. At the transition point the distribution should be 
independent of the system size. In Fig. [5] we plot the 
behavior of 7 as function of W for several sizes of a scale 
free graph. Indeed, 7 decreases with system size for small 
values of W while it increases with size for large values 
of M^. 

All curves should cross around a particular value of 
disorder signifying the critical disorder. From finite size 
scaling arguments ^30j one expects that 7 around the crit- 
ical disorder will depend on the the disorder and network 
size, L, in the following way: 



Pp{s) 



(5) 



^{W,L)^^{W„L) + C 



W 



1 



(7) 



4 



f~ 0, 




FIG. 2: 7 as function of W for different SF graphs sizes (A = 4 
and m = 2). The typical behavior for finite size transition is 
seen, where a crossing in the size dependence of 7 between 
the metallic (small values of W) and localize (large value of 
W) regime is seen. 
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FIG. 3: The scaling of 7 according to Eq. for different SF, 
A = 4, m = 2, networks sizes. Two branches, corresponding 
to the metallic and localized regimes, appear. 



where C is a constant. This relation enables us to extract 
both the critical disorder Wc and the critical index v. 
Scaling of the numerical data according to Eq. {T)) yields 
two branches corresponding to the metallic and localized 
regimes, that are clearly seen in Fig. |3 The estimated 
values of v and Wc (see Table I) are extracted by fitting 
the branches to a 4th order polynomial. 



Network 


(k) 


I 


Wc 


V 


Scale-Free , A = 4 , m = 2 


2.97 


12.46 


15.7±0.9 


0.55±0.11 


Random- Regular (RR) 


3 


11.8 


11.9±0.26 


0.66±0.08 


Erdos-Renyi 


3 


9.45 


20.5±0.23 


0.68±0.08 


Cayley-Tree 


3 


10 


11.44±0.06 


0.51±0.045 


Cayley-Tree with loops 


3.05 


10 


12.4±0.1 


0.54±0.075 


RR "double peak" 


3.1 


10.28 


14.1±0.3 


0.85±0.41 



TABLE I: Networks showing the localization transition. The 
value of / is for N = 1000. 



III. RESULTS AND DISCUSSION 

The calculations for all networks mentioned above 
are performed for M different realizations, where 
M = 1000,400,200, ...,50 for the corresponding num- 
ber of nodes: N = 200, 500, 1000, 4000. Ex- 
cept for the Cayley-tree networks for which M — 
4000, 2000, 1000,..., 125, 64 for the corresponding tree 
sizes: TV = 63, 127, 255, 2047, 4095 or L = 
6, 7, 8, ... , 11, 12 (where L is the number of "generations" 
of the tree). Another exception is for Erdos-Renyi graphs 
in which (/c) is between 3 and 3.5. The low connectivity 
of the graphs, results in one main cluster and relatively 
large number of not-connected nodes (about 5%). Thus, 
the calculations are made only for the largest cluster of 
each realization, since a procedure that considers all the 
nodes is skewed by the eigenvalues of small disconnected 
clusters [sj. A clear localization transition is observed 
for a group of graphs which are all characterized by an 
average degree (/c) smaller than 3.1, and an averaged last 
occupied shell I (for N — 1000 sites) larger or equal to 
9.45. The results are summarized in Table I. 

The results for all the graphs (including those which 
show no clear signs of transition) can be scaled according 
to their average degree (fc). The higher the value of (A;) is 
- the higher is the value of W needed in order to obtain a 
specific value of 7. Thus, the higher the average degree, 
the more metallic the system is, which makes sense. A 
cross section at 7 = 0.6 of all curves is shown in Fig0| 
as a function of (fc). The (fc) of the networks studied in 
Fig. ^ as well as the averaged last occupied shell I for 
N = 1000 sites are presented in Table II. 

The following observations can be gleaned out of the 
data for the different networks: 

(1) For all the networks that show a metal-insulator tran- 
sition, V is of order 1 /2 except for the Random- Regular 
" double-peak" network which is the one with the highest 
value of connectivity that still shows a clear transition. 
A critical index of = 0.5 is expected for a system of 
infinite dimensionality. At (fc) = 3.1 the value of v is 
significantly higher, but so is the estimate of the error 
bar. On the other hand for the Erdos-Renyi graph with 
(/c) =3.1 no clear transition is observed. 

(2) All networks with connectivity above 3.1 do not show 
clear signs of a metal-insulator transition. Nevertheless, 
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Network 


(k) 


I (for iV = 1000) 


Scale-Free , A = 4 , m = 2 


2.97 


12.46 


Random- Regular 


3 


11.8 




Random- Regular "double peak' 








iP 


= 0.95 -^k = 3 ,p = 0.05 k = 


= 5) 


3.1 


10.28 


Cayley-Tree 


3 


10 


Cayley-Tree with loops 


3.05 


10 


Erdos-Renyi 


3 


9.45 


Scale-Free , A = 3.5 , m = 2 


3.28 


9.36 


Erdos-Renyi 


3.1 


9.31 


Erdos-Renyi 


3.2 


9.03 




Random- Regular "double peak' 










) = 0.8 ^ fe = 3 , p = 0.2 fc = 


5) 


3.4 


9 


Erdos-Renyi 


3.5 


8.33 




Random- Regular "double peak' 








{p-- 


= 0.95 ^ fc = 3 , p = 0.05 ^ = 


10) 


3.35 


7.99 


Erdos-Renyi 


4 


7.51 


Scale-Free , A = 5 , m = 3 


4 


7.37 


Scale-Free , A = 4 , m = 3 


4.5 


6.05 




Random- Regular "double peak' 








(P 


= 0.8 ^ fc = 3 , p = 0.2 ^ fc = 


10) 


4.4 


6.13 


Erdos-Renyi 


5 


6.31 


Erdos-Renyi 


7.5 


5.02 


Erdos-Renyi 


10 


4.1 



TABLE II: The average connectivity (k) of all the networks 
considered in this study, as well as the averaged last occupied 
shell I ior N ^ 1000 sites. 
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one should be rather careful in interpreting this observa- 
tion since, as is clear from Table II, larger values of (k) 
lead to smaller size, /, of the network for the same number 
of nodes. Moreover, from the two networks which have 
the same (fc) =3.1, only the one with the higher value 
of I shows clear signs of the metal insulator transition. 
Thus, the absence of transition may be an artifact of the 
small size of networks with high average connectivity. 

(3) The critical disorder Wc fluctuates in the range of 
12-20 (Table I). Due to the small range of (k) (2.97-3.1), 
it is hard to determine any relation between k and Wc- 

(4) On the other hand, there is a clear relation between 
the amount of disorder needed in order to reach a par- 
ticular value of 7 (i.e., the value of W needed to reach 
a certain degree of localization) and (k) . As can be seen 
in Fig. ^ a linear dependence Wi^j = 0.6) cx {k) is ob- 
served. 



Thus, the gross features of the Anderson metal- 
insulator transition are similar for a wide range of dif- 
ferent networks. The critical index for all the networks 
studied here are within the range expected for a system 
of infinite dimensionality, and the connectivity influences 
the degree of localization. Thus, complex networks are an 
example of a topology for which the critical index follows 
the mean field prediction, but the critical disorder can be 
tuned by the connectivity. On the other hand, the fact 
that networks with high connectivity are very compact 
raises the problem of identifying the transition point. It 
is hard to extend the usual finite-size scaling method to 
networks with high connectivity since the number of sites 
grows very rapidly with size, while for small network sizes 
the crossover behavior of the 7 curves is very noisy. This 
results in an inability to clearly identify the Anderson 
transition, although it can not be ruled out the possibility 
that there is a critical connectivity for complex networks 
above which no transition exist. 



20 
10 



6 

<k> 
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FIG. 4: The 7 curves of the checked networks can be ordered 
by their (k). The values of their W{'y — 0.6) are presented as 
a function of (k). One can notice the increasing W with (k). 
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